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Abstract 

Many astrophysical data show that the expansion of our universe is accelerating. 
In this paper, we study the model of phantom with O(N) symmetry in background 
of loop quantum cosmology(LQC). We investigate the phase-space stability of the 
corresponding autonomous system and find no stable node but only 2 saddle points 
in the field of real numbers. The dynamics is similar to the single-field phantom 
model in LQC[1|. The effect of O(iV) symmetry just influence the detail of the 
universe's evolution. This is a sharp contrast with the result in general relativity, in 
which the dynamics of scalar fields models with O(N) symmetry are quite different 
from the single-field models [2J,[3J,[4J. 
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1 Introduction 



Recently, many astrophysical data, such as WMAPj5], type la supernovae[6] and large 
scale galaxy surveys [7], show that our universe is almost flat, and currently undergoing 
a period of accelerating expansion which is driven by a yet unknown dark energy. About 
70% of the energy density in our universe is dark energy nowadays. There are many 
candidates of dark energy such as the cosmological constant [8], quintessence [9] |10j CHI 
IT2| IT3~t |2], phantomJHl [151 HHJ E], quintom[T7J H] etc. The equation-of-state parameter 
w = p/p plays an important role in these models. The parameter of quintessence and 
phantom lies in the range — l/3>u»>— 1 and w < — 1 respectively. And for quintom w 
could cross -1. 

Li et al. generalized the quintessence and phantom models to fields possess O(iV) 
symmetry[2] [3]. Setare et al. extended the generalization to the quintom model[4j. In 
the above papers they showed that the behaviors of the dynamic system with O(iV) 
symmetry are quite different from the single-field models. 

Many works on dark energy have been done in the framework of Einstein's classical 
general relativity(GR). However, most physicists believe that the gravity should be quan- 
tized in a ultimate theory. Thus, we should consider the quantum effects in the evolution 
of our universe. Loop quantum gravity(LQG) is a outstanding candidate of quantum 
gravity theory. LQG is a background independent and non-perturbative theory. At 
the quantum level, the classical spacetime continuum is replaced by a discrete quantum 
geometry and the operators corresponding to geometrical quantities have discrete eigen- 
values. Many topics on cosmology are discussed within the framework of LQG in various 
literature in which it is known as Loop Quantum Cosmology (LQC) [H [18l [19l EQl [21]. Re- 
cently, Samart et aLpQ have studied the phantom dynamics in LQC. The work of quintom 
and hessence model in loop quantum cosmology have been done by Wei et al.\V3\. They 
found the behaviors is different from the ones in standard FRW, such as the avoidance 
of future singularities. 

In this paper, we investigate the evolution of universe dominated by phantom with 
O(iV) symmetry. The phase-space stability of the corresponding autonomous system has 
been discussed. We find there's no stable attractor but only 2 saddle points, this result 
is similar to the single-field phantom in LQC. In order to research the evolution of our 
model, we numerically study the system in exponential potential, and draw 2 graph to 
display the evolution of Hubble parameter H and energy density p to cosmic time t. 
The figures reveal that in LQC even for a universe dominated by phantom, the Hubble 
parameter H may turn negative and oscillate. Thus, the universe will enter a oscillatory 
regime. The relation between p and H is also shown in a graph. 
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2 Loop quantum cosmology 



Loop quantum cosmology(for review, see Ref . [221 [23] ) can describe the homogeneous 
and isotropic spacetimes by effective modified Friedmann equation for flat universe. The 
quantum effects could be reflected by adding a correction term p/p c into the eqution. 
In Einstein's general relativity, the standard Friedmann equation which describes a flat 
universe is: 

H 2 = yP (2.1) 
Instead of Eq. ( 12. ip . the modified Friedmann equation for a flat universe in LQC is given 

by (snug 

H 2 = yP(l - ±) (2.2) 
3 p c 

As in GR, H = a/a is Hubble parameter, p is the total energy density, and a dot denotes 
the derivative with respect to cosmic time t. p c is critical loop quantum density: 

Pc = 16ttW^' (2 - 3) 

Where 7 is the dimensionless Barbero-Immirzi parameter and is suggested that 7 ~ 
0.2375 by the black hole thermodynamics in LQGj23]. The modified Friedmann equation 
provides an effective description for LQC which very well approximates the underlying 
discrete quantum dynamics. Differentiating Eq. (12. 2p and use the conservation law 

p + 3H (p + p) =0 (2.4) 

We obtain the effective modified Raychaudhuri equation 

H=~{p + p)(l-2^) (2.5) 

* Pc 

where p is the total pressure. In LQC the accelerating condition is 

" H + H 2 >0 (2.6) 



a 



The dynamics of our universe could be analyzed through Eq. (12.21) . (12 .4p . (12. 5p and the 
equation of the fields which dominate the universe. 
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3 O(N) Phantom 



We consider the flat Robert son- Walker metric 

d s 2 = dt 2 - a 2 (t) (dx 2 + dy 2 + dz 2 ) (3.7) 
The Lagrangian density for the Phantom with O(N) symmetry is given by[3] 

U = ~g» v {d^ a ){d v ® a ) - (3.8) 

where <3> a is the component of the scalar field, a = 1, 2, • • ■ , N. In order to impose the 
O(N) symmetries, following[3j, we write it in the form 

$ x = R(t) cosy^t) 
$ 2 = R(t) siny?i(t) cosy2 2 (t) 
$ 3 = R(t) sin (p!(t) sin (p 2 (t) cos ip 3 (t) (3.9) 



$ N 1 = R(t) sin(pi(t) ■ ■ ■ sintpN~2(t) cosipN-i(t) 
& N = R(t) sin (fi(t) ■ ■ ■ sin(pN~2(t) sin (fN-i(t) 

Thus, we have |$ a | = R and assume that the potential V(|$ a |) depends only on R. 
The radial equation of scalar fields is 

The "angular components" contribute a effective term Jrp- to the system's dynamics. 
The energy density p and the pressure p of the O(N) phantom are given as 

P = ~\ {k2 + ^ ) + V{R) (3 ' n) 

and 

Hence, the equation-of-state for the O(N) phantom is 

w = -p-. — 3.13 

Different from GR, the accelerating condition for universe in LQC is not w < —1/3 but 
Eq.fl22D. 
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4 EVOLUTION AND STABILITY OF THE MODEL 



In this section we discuss the dynamics of universe dominated by phantom with O(N) 
symmetry in LQC. In order to get a possible evolution of the model, we choose exponential 
potential 

V = V e- XKR , (4.14) 

The evolution of universe is governed by Eq. (12 . 2 1) . (12 . 5 j) . ( 13TT0|) . To investigate the stability 
of the model, we introduce the following dimensionless variables 



x = J*- V =^/XK Z =^JL (415) 
m = —, f = — , iV = lna (4.16) 

Using these variables, Eq. ( 12. 2p . (12. 5p . (13. 10p can be written as the following autonomous 
system: 



-3x + VQ£z 2 - ^Ay 2 - (3a; 3 + 3xz 2 )(l - 2m) 



dx 
dN 

d V 1 3 ... x .2 ,,. .2 



3z- v / 6x^-(3x 2 ^ + 3z 3 )(l-2m) (4.17) 



rfA^ V 2 

dN 

dm —x 2 — y 2 — z 2 . 

_ = _ 3m(i + _____ 
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In the following, it is straightforward to analyze the critical points as well as their 
stability. We can get critical points by setting the right hand of the equation (l4.17p to 
zero. As presented in table [U in the field of real numbers there are only 2 real critical 
points: A and B. We expand the variables around the critical point A and B in the form 
x = x c + Sx, y = y c + 6y, z = z c + 5z,m = m c + 6m, £ = £ c + 5£, where Sx, Sy, Sz, 5m, 5£ 
are perturbations of the variables near the critical points and we consider them forming a 
column vector denoted as U. Substituting the above expression into (I4.17P and neglecting 
higher order terms in the perturbations, we can obtain the equations for the perturbations 
up to first order as: 

U' = M • U, (4.18) 

where the prime denotes differentiation with respect to N. M is a 5 x 5 matrix formed 
by the coefficients of the perturbation equations (I4.18p . The matrix for the critical point 
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Critical point 


x c 


Vc 




m c 




A 


X 


V6+A 2 
J% 











B 


X 


V6+A 2 












Table 1: The real critical points of the autonomous system f)4.17p . 



Label 


Eigenvalues 


Stability 


w 


Acceleration 


A 


- X \ |(-6-A 2 ), |(-6-A 2 ) , A 2 ,0 


saddle point 


I 


Yes 


B 


-A 2 , ±(-6 -A 2 ), ±(-6 -A 2 ), A 2 , 


saddle point 


-1 - ^ 

3 


Yes 



Table 2: The eigenvalues of the matrix M and the properties of critical points in the 
autonomous system (I4.17p . 



(x c , y c , z c , m c , £ c ) is given by 



/ dx' dx' dx' dx' dx' 

dx dy dz dm d£ 

dy' dy' dy 1 dy' dy' 

dx dy dz dm 9§ 

dz' dz' dz' dz' dz' 

dx dy dz dm d£ 

dm' dm' dm' dm' dm' 

dx dy dz dm d£ 

d? d£' dj' dj' 

\ dx dy dz dm d£ / ^ x=Xciy=ycjZ=Zc , m=mc ^ = ^ 



The eigenvalue of (I4.19P determine the type and stability of the critical points. We 
present them in table [2j As shown, both critical A and B are saddle points, it follows 
immediately that the phase trajectory is very sensitive to initial conditions given to the 
system. This result is similar to the single- field phantom in LQC[lj. It is interesting to 
contrast the influence of 0(N) symmetry in GR and LQC. Several worksj2]|3]Il] shows 
that the dynamics of quintessence, phantom and quintom with 0(N) symmetry in GR 
is quite different from the single-field model. 



(4.19) 



In the following, we study the possible evolution of our model numerically. In order 
to do this, we set the value of some parameters: A = —0.1, Q = 1, k = V8n, p c = 
1.5, Vo = 1.2, and draw 3 graphs: FigJH Figj2j Figj3j FigJT] shows the evolution of 
Hubble parameter H with respect to cosmic time t. From Fig{TJ we found H may take 
maximum value while p = p c /2, and then decline even to negative interval. After H taking 
minimum value, it would bounce and then oscillate. Thus, the accelerating expansion of 
our universe would end in finite time and undergoes contraction from halting(if = 0). 
Then the universe enters an oscillating stage. Figj2] shows the evolution of cosmic total 
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energy density p, it oscillate while t goes by. In fact, from FigJ3]we can get some insight 
into the evolution of the model. FigJ3]reveals the relation between H and p. As time goes 
by, the trajectory makes a closed region, and the oscillation of H and p are restricted in 
a certain interval. 




Figure 1: The evolution of H with respect to t for A = —0.1, Q — 1, k — \/8tt, p c = 1.5, 
V = 1.2. 




8 9 10 11 12 13 14 t 

Figure 2: The evolution of cosmic total energy density p with respect to t for A = —0.1, 
n = l,K = v 7 ^, p c = 1.5, Vb = 1.2. 



5 Conclusions 



In this paper, we investigate the dynamics of phantom model with O(N) symmetry 
in background of loop quantum cosmology. We discuss its stability by analyzing the 
autonomous system Eq. fl4.17p . It shows that there is no attractor in the system, but 
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. 



-0 . 



-1 . 

Figure 3: The relation between H and p while t goes from to 14. Set values are 
A = -0.1, n = l,K= V8^, p c = 1.5, V = 1.2. 

only 2 saddle points. This result is similar to the single-field model in LQC([T]). The 
effect of O(N) symmetry just influence the detail of the universe's evolution. It is a sharp 
contrast with the result in GR, in which the dynamics of scalar fields models with O(N) 
symmetry are quite different from the single-field models. Finally, to show some insight 
into our model, we draw 3 graphs by Numerical analysis. The possible oscillation of H 
and p are revealed in Fig{TJ and Figj2] Figj3] shows that as time goes by, H and p are 
restricted in a certain interval, and the trajectory makes a closed region. 
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